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Abstract 

This article characterizes the singularities of very weak solutions of 3D sta- 
tionary Navier-Stokes equations in a punctured ball which are sufficiently small 
in weak L 3 . 

Keywords: stationary Navier-Stokes equations, point singularity, very weak 
solution, Landau solution. 



1 Introduction 

We consider point singularities of very weak solutions of the 3D stationary Navier- 
Stokes equations in a finite region Q in R 3 . The Navier-Stokes equations for the 
velocity u : Q — > R 3 and pressure p : Q — > R with external force / : Q — > R 3 are 

- Au + (u ■ V)u + Vp = /, divu = 0, (i G fl). (1.1) 

A very weak solution is a vector function u in Lf oc (Q) which satisfies (11.11) in distri- 
bution sense: 

-u ■ A<p + u jUi d j(Pi = (f, <p), V(/j E CZ(ty, (1.2) 



and J u ■ Vh = for any h G C^°(Q). Here the force / is allowed to be a distribution 
and 

C~(f2) = {<p e C c °°(ft,R 3 ) : div V = 0}. (1.3) 
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In this definition the pressure is not needed. Denote Br = {x G M. 3 : \x\ < R} and 
B% = K 3 \ B R for R > 0. 

We are concerned with the behavior of very weak solutions which solve ( 11.11) in 
the punctured ball B 2 \{0} with zero force, i.e., / = 0. There are a lot of studies 
on this problem [5], [TU [T5], HI HO] . A typical result is to show that, under some 
conditions, the solution is a very weak solution across the origin without singular 
forcing supported at the origin (removable singularity), and is regular, i.e., locally 
bounded, under possibly more assumptions (regularity). Dyer-Edmunds [5] proved 
removable singularity and regularity assuming both u,p G L 3+£ (B2) for some e > 0. 
Shapiro [T4"l [15] proved removable singularity and regularity assuming u G L 3+e (B2) 
for some e > and u(x) = odx)" 1 ) as x — > 0, without assumption on p. Choe and 
Kim jl] proved removable singularity assuming u G L 3 (B 2 ) or u{x) = o(|x| _1 ) as 
x — > 0, and regularity assuming it G L 3+e (5 2 ) for some e > 0. Kim and Kozono 
[TO] recently proved removable singularity under the same assumptions as jl], and 
regularity assuming u G L 3 (B 2 ) or m is small in weak L 3 . As mentioned in [TU], their 
result is optimal in the sense that if their assumption is replaced by 

|u(x)| < C.\x\~ x (1.4) 

for < | a; | < 2, then the singularity is not removable in general, due to the existence 
of Landau solutions, which is the family of explicit singular solutions calculated by 
L. D. Landau in 1944 [8], and can be found in standard textbooks, see e.g., j£l p. 82] 
or P p. 206]. 

The purpose of this article is to characterize the singularity and to identify the 
leading order behavior of very weak solutions satisfying the threshold assumption 
(11. 4p when the constant is sufficiently small. We show that it is given by Landau 
solutions. 

We now recall Landau solutions in order to state our main theorems. Landau 
solutions can be parametrized by vectors b G M. 3 in the following way: For each 
feel 3 there exists a unique (— l)-homogeneous solution U b of (II. ip together with an 
associated pressure P b which is (— 2)-homogeneous, such that U\ P b are smooth in 
M 3 \{0} and they solve 

- Ait + (u ■ V)it + Vp = bS, div u = 0, (1.5) 

in M 3 in the sense of distributions, where 5 denotes the Dirac 5 function. When b = 
(0, 0, (3) with f3 > 0, they have the following explicit formulas in spherical coordinates 
r, 9, <p with x = (r sin 9 cos0, r sin 9 sin <p, rcos9): 

2/ A 2 -l \ 2 sing -4(Aco S fl-l) 

"rV(^-cos#) 2 ) 6r r(A-cos9) ed) r 2 (A-cos^) 2 1 } 
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where e r = ~ and eg = (cos 9 cos (f>, cos 6* sin 0, — sin^). The parameters j3 > and 
A G (1, oo] are related by the formula 

^16,(^ + -^lo g _ + ^- TI j. (1.7) 

The formulas for general b can be obtained from rotation. One checks directly that 
rf/ b || Loo is monotone in \b\ and || r ^ fe || ioo ~~ ¥ (or oo) as \b\ — > (or oo). Recently 
Sverak [16] proved that Landau solutions are the only solutions of (II. ip in IR 3 \{0} 
which are smooth and (— l)-homogeneous in M 3 \{0}, without assuming axisymmetry. 
See also [TSl [21 E] for related results. 

If u,p is a solution of (1.1), we will denote by 

T i:j (u,p) = p5 i:j + itjitj - - djUi (1.8) 

the momentum flux density tensor in the fluid, which plays an important role to 
determine the equation for (u,p) at 0. Our main result is the following. 

Theorem 1.1 For any q G (1,3), there is a small C* = C*(q) > such that, if u is 
a very weak solution of (11.11) with zero force in _B 2 \{0} satisfying (11.41) in I?2\{0}, 
then there is a scalar function p satisfying \p{x)\ < C\x\~ 2 , unique up to a constant, 
so that (u,p) satisfies (11.51) in B 2 with b{ = r, =1 T i j(u,p)rij(x) ) and 

\\ u ~ ub \\w^(BA + SU P kl 3/9_1 |(^ " U b )(x)\ < CC*, (1.9) 

where the constant C is independent of q and u. 

The exponent q can be regarded as the degree of the approximation of u by U b . 
The closer q gets to 3, the less singular u — U b is. But in our theorem, C*(q) shrinks 
to zero as q — > 3_. Ideally, one would like to prove that u — U b G L°°. However, it 
seems quite subtle in view of the following model equation for a scalar function, 

-Av + cv = 0, c = Av/v. (1.10) 

If we choose v = log |xj , then c(x) G L 3 ^ 2 and lim| x |^ l a; | 2 | c ( a; )l = 0, but v G" L°°. 
In equation (13.21) for the difference w = u — U b , there is a term (w ■ V)U b which has 
similar behavior as cv above. 

In fact, we have the following stronger result. Denote by U wk the weak U spaces. 
We claim the same conclusion as in Theorem 11.11 assuming only a small L z mk bound 
of u but not the pointwise bound (11.41) . 
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Theorem 1.2 There is a small 5* > such that, if u is a very weak solution of 
(11. ip with zero force in Q = B 21 \{0} satisfying \\u\\ L 3 (n) =: e < £*, then u satisfies 
\u[x)\ < Ci^lxl" 1 in B 2 \{0} for some C\. Thus the conclusion of Theorem \l.l\ holds 
ifC x e<C*(q). 

Our results are closely related to the regularity problem of very weak solutions, 
which could be considered when u is only assumed to be in Lf oc . In fact, the problem 
with the assumption u being large in L 3 ^ already exhibits a great difficulty Recall 
the scaling property of (11. 1ft : If iu,p) is a solution of (ll.ip . then so is 

(u x ,p x )(x) = (A«(Ax),A 2 p(Ax)), (A > 0). (1.11) 

The known methods are primarily perturbation arguments. Since L 3 fe -quasi-norm 
is invariant under the above scaling and does not become smaller when restricted to 
smaller regions, one would need to exploit the structure of the Navier-Stokes equations 
in order to get a positive answer. Compare the recent result [3] on axisymmetric 
solutions of nonstationary Navier-Stokes equations, which also considers a borderline 
case under the natural scaling. 

This work is inspired by Korolev-Sverak [11] in which they study the asymptotic 
as |x| — > oo of solutions of (II. ip satisfying ( j 1.4ft in M 3 \i?i. They show that the 
leading behavior is also given by Landau solutions if is sufficiently small. Our 
theorem can be considered as a dual version of their result. However, their proof is 
based on the unique existence in M 3 of the equation for v = ipiu — U b ) + £ where ip is 
a cut-off function supported near infinity and ( is a suitable function chosen to make 
div-u = 0. If one tries the same approach for our problem, since one needs to remove 
the origin as well as the region |x| > 2 while extending u — U b , one needs to choose a 
sequence ifk with the supports of 1 — ifk shrinking to the origin, which produce very 
singular force terms near the origin. 

Instead, we first prove Lemma 12.31 which gives the equation for (u, p) near the 
origin. Since the equation for u is same as U b near the origin for b = b(u), the <5- 
functions at the origin cancel in the equation for their difference. We then apply the 
approach of Kim-Kozono [10] to the difference equation, and prove its unique existence 
in Wl ,r {B 2 ) for 3/2 < r < 3 and uniqueness in W 7 1 ' r nL 3 ofc (_B 2 ) for 1 < r < 3/2, which 
improves the regularity of the original difference. Above Wq' v (B 2 ) is the closure of 
C c °°(5 2 ) in the W 1 ' r (B 2 )-norm. 

As an application, we give the following corollary. Recall u\ for A > is defined in 
(II. lip . A solution u on £> 2 \{0} is called discretely s elf- similar if there is a Ai e (0, 1) 
so that u\ 1 = u. Such a solution is completely determined by its values in the annulus 
Bi\B\ 1 since u(\\x) = Xi k u(x). They contain minus-one homogeneous solutions as 
a special subclass. 
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Corollary 1.3 If u satisfies the assumptions of Theorem \1.1\ and furthermore u is 
discretely self-similar in B 2 \{0}, then u = U b . 

This corollary also follows from [TT] (with domain M 3 \i?i and Ai > 1). In the 
case of small C*, this corollary extends the result of Sverak [16] on minus-one homo- 
geneous solutions. The classification of discretely self-similar solutions with large C* 
is unknown. 

As another application, we consider a conjecture by Sverak [TBI §5]: 

Conjecture 1.4 If u is a solution of the stationary Navier-Stokes equations (11.11) 
with zero force in M 3 \{0} satisfying ( 11.41) with some C* > 0. Then u is a Landau 
solution. 

We give a partial answer for this problem. 



Corollary 1.5 Conjecture\l.J\ is true, provided the constant C* is sufficiently small. 



The above corollary can be also shown to be true by either our main theorem or 
the result of Korolev-Sverak [TT], see section 13741 The corresponding conjecture for 
large C* is related to the regularity problem of evolutionary Navier-Stokes equations 
via the usual blow-up procedures. 



2 Preliminaries 

In this section we collect some lemmas for the proof of Theorem ll.il The first lemma 
recalls Holder and Sobolev type inequalities in Lorentz spaces. We denote the Lorentz 
spaces by L p ' q (1 < p < oo, 1 < q < oo). Note L 3 ^ = L 3 '°°. 

Lemma 2.1 Let B = B 2 C W 1 , n > 2. 

i) Let 1 < pi,p2 < oo with 1/p := l/p\ + l/p 2 < 1 and let 1 < ri,r 2 < oo. For 
f E L Pl ' ri and g E L P2 ' r2 , we have 

\\fg\\LP,r(B) < C || /IU?i. -1(B) || #11^2. '-2(B) for r := min{r 1 ,r 2 }, (2.1) 

where C = C(p 1 , ri,p 2) r 2 ). 

ii) Let 1 < r < n. For f G W 1,r (B), we have 

\\f\\ L ^r [B) <C\\f\\ W l,r {B) , (2.2) 

where C = C(n,r). 
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Part (i) of Lemma [2.11 was proved in [13]. Part (ii) was proved in [13] for M. n and 
in [12l [10] for bounded domains. 

By this lemma, when n = 3 and 1 < r < 3, we have 

\\fg\\ L r {B) < C \\f\\ L s wk \\g\\ L &, r < C r \\f\\ Llk{B) \\g\\ w i,r {B) . (2.3) 

This estimate first appeared in [10] and plays an important role for our application. 

The next lemma is on interior estimates for Stokes system with no assumption on 
the pressure. 

Lemma 2.2 Assume v £ L 1 is a distribution solution of the Stokes system 

- Avi + = djfij, div v = in B 2R (2.4) 

and f £ U for some r £ (1, oo). Then v £ W^ r c and, for some constant C r indepen- 
dent of v and R, 

\\Vv\\ Lr{BR) < C r \\f\\ Lr{B2R) + C r R~^ r \\v\\ Ll(B2R) . (2.5) 

This lemma is (T7J, Theorem 2.2. Although the statement in [T7] assumes v £ 
W l( Jc , its proof only requires v £ L 1 . This lemma can be also considered as [HI Lemma 
A. 2] restricted to time- independent functions. 

The following lemma shows the first part of Theorem ll.il except (11.91) . In partic- 
ular, it shows that (u,p) solves (11. 5p . 

Lemma 2.3 If u is a very weak solution of (11.11) with zero force in B 2 \{0} satisfying 
(ll.4p in i? 2 \{0} (with C* allowed to be large), there is a scalar function p satisfying 
\p(x)\ < C\x\~ 2 , unique up to a constant, such that (u,p) satisfies (II. 5p in B 2 with 
h = J^ =1 Tij(u,p)nj(x) . Moreover, u,p are smooth in B 2 \{0}. 

Proof. For each R £ (0, 1/2], u is a very weak solution in B 2 — B B in L°°. Lemma 
12.21 shows u is a weak solution in W^. The usual theory shows that u is smooth and 
there is a scalar function p B , unique up to a constant, so that (u,p B ) solves (11.11) in 
B 2 — Br, see e.g. [7]. By the scaling argument in Sverak-Tsai [17] using Lemma |2~2| 
we have for x £ B 3R — B 2R , 

\V k u(x)\<^^ for k = 1,2,..., (2.6) 

where = C/^C*) are independent of R £ (0, 1/2] and its dependence on can be 
dropped if C* £ (0, 1). Varying R, §ZM) is valid for x £ B 3/2 \{0}. For < R < R', 
by uniqueness of p' R , the difference Pr\ B2 - B , — Pr> is a constant. Thus we can fix 
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the constant by requiring pa = py 2 in B 2 \Bi/ 2 , and define p(x) = Pr(x) for any 
x G B 2 \{0} with R = \x\/2. By the equation, |Vp(x)| < CC*\x\~ 3 . Integrating from 
\x\ = 1 we get \p{x)\ < CC* \x\~ 2 . In particular 

|7ii(«,p)(a0| < CC*\x\~ 2 for x G B 3/2 \{0}. (2.7) 

Denote NS(u) = —Au + (u ■ V)n + Vp. We have NS(u)i = djTij(u) in the sense 
of distributions. Thus, by divergence theorem and NS(u) = in B 2 \{0}, 



h= T ij (u,p)n j (x) = / 7ij(«,p)nj(a;) (2.8) 

J|x|=l J|x|=_R 

for any i? G (0, 2). Let be any test function in C%°(Bi). For small s > 0, 

(NS(u)i,<j>) = - [ T l3 (u)d 3 <P 



djTij(u)(f)+ / Tij(u)(j)nj - / Tij(u)(jmj— / Tij{u)dj(j). 

Bi\B e JdB E JdBx J B e 

In the last line, the first integral is zero since NS(u) = and the third integral is zero 
since = 0. By the pointwise estimate ( 12. 71) . the last integral is bounded by Ce 3 ~ 2 . 
On the other hand, by ( 12.81) . 

Tij(u)<jmj -> M>(0) as e -> 0. (2.9) 

Thus (w,p) solves (11.51) and we have proved the lemma. □ 

It follows from the proof that \b\ < CC* for C* < 1. With this lemma, we have 
completely proved Theorem 11.11 in the case q < 3/2. In the case 3/2<g<3, it 
remains to prove ( 11.91) . 

3 Proof of main theorem 

In this section, we present the proof of Theorem 11.11 We first prove that solutions 
belong to W 1,q . We next apply this result to obtain the pointwise estimate. For what 
follows, denote 

w = u-U, U = U\ (3.1) 
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where U b is the Landau solution with b given by (12.81) . By Lemma 12.31 there is a 
function p such that (w,p) satisfies in B 2 that 

-Aw + U ■ Vw + w ■ V(U + w) + Vp = 0, div w = 0, 

i / m ^ i -/ m - (3-2) 

Note that the ^-functions at the origin cancel. 
3.1 W 1,q regularity 

In this subsection we will show w G W 1,q (Bi). Fix a cut off function ip with tp = 1 in 
5 9 /8 and </? = in B^ g . We localize w by introducing 

u = <y9w + £ (3.3) 

where £ is a solution of the problem div£ = — • w. By Galdi 0, Ch.3] Theorem 
3.1, there exists such a ( satisfying 

suppC C £ 3 / 2 \£i, HVClUioo < C||V^-w|| L ioo < CO, (3.4) 

The vector w is supported in -B3/2, satisfies v G jy 1 ^ D for r < 3/2 by (II .41) . (12.61) 
and (13.41) . and 

- Av + U ■ Vv + v ■ V(U + v) + Vvr = /, div v = 0, (3.5) 
where ir — <pp, and 

/ = -2(V(f ■ V)w - (Atp)w + (U ■ V(f)w + (if 2 - ip)w -Vw + (w- V(f)w 



+pV<p - AC + (U ■ V)C + C • V(f/ + + C) + ■ VC 
is supported in the annulus B 3 / 2 \Bi. One verifies directly that, for some C\, 



(3.6) 



sup \\f\\ w -^ m < CiC*. (3.7) 

l<r<100 

Our proof is based on the following lemmas. 

Lemma 3.1 (Unique existence) For any 3/2 < r < 3 ; for sufficiently small C* = 
C„(r) > 0, there is a unique solution v of (13. 5p and (13.71) in the set 

V = {v G W 1 ' r (£ 2 ), H y := IMI^r^j < C 2 a} (3.8) 

/or some C 2 > independent of r G [3/2,3). 



8 



Lemma 3.2 (Uniqueness) Let 1 < r < 3/2. If both v\ and v 2 are solutions of 
(13. 5p and (13.71) in Wq' t fl and C* + \\vi\\ L 3 + IKH^ sufficiently small, then 

wk wk 

Vi = v 2 . 

Assuming the above lemmas, we get W 1,q regularity as follows. First we have a 
solution v of (13.51) in WQ ,q (B 2 ) by Lemma [3.11 On the other hand, both v = ipw + ( 
and v are small solutions of ( 13. 5\\ in Wq 1 " fl Ll, k (B 2 ) for r = 5/4, and thus w = v by 
Lemma E2 Thus v G w£' q {B 2 ) and w G W^-Bi). 

Proof of Lemma 13.11 Consider the following mapping $: For each w G V, let 
w = $-y be the unique solution in W Q ' r (B 2 ) of the Stokes system 

- Au + V7f = /- V- (U®v + v® (U + v)), divI7 = 0. (3.9) 

By estimates for the Stokes system, see Galdi [6j Ch.4] Theorem 6.1, in particular 
(6.9), for 1 < r < oo, we have 

\\v\\ w t,r m < N r \\f\\ w -i,r + N r \\V-(U®v + v®(U + v))\\ w -i,r (3.10) 

for some constant N r > which is uniformly bounded for r in any compact regions 
of (1, oo). By (13. 7p and Lemma [2TT1 in particular (I2.3p . for 1 < r < 3, 

IMIwo 1 '^) ^ N r CxC* + N r \\U 0v + v0(U + v)\\ L r 

< N r da + N r C r (\\U\\ L S wk + \\v\\ Llk )\\v\\ V . (3,U) 

We now choose C 2 = 2(Ci + 1) sup 3 / 2<r<3 N r . Since V C Lf uk if r > 3/2, we get 
v = $w G V if C* is sufficiently small. 

We next consider the difference estimate. Let v±, v 2 G V, Vi — and w 2 = $u 2 . 
Then 

ll^'i - ^v 2 \\ w i,r < CC r {\\U\\ L a wk + IkilU^ + IHLaJHu! - V 2 \\ W l,r. (3.12) 

Taking C* sufficiently small for 3/2 < r < 3, we get \\&Vi — &v 2 \\ v < \ \v\ — v 2 \\ v , 
which shows that $ is a contraction mapping in V and thus has a unique fixed point. 
We have proved the unique existence of the solution for (I3.5l) - (l3.7p in V. □ 

Remark. Since the constant C r from Lemma \2. II (ii) blows up as r — > 3_, our C* 
shrinks to zero as r — > 3_. 

Proof of Lemma 13.21 By the difference estimate (13.121) . we have 

IK -v 2 \\ w i,r < C(\\U\\ L Z + \\Vl\\ L 3 + 1 1 ^2 1 1 i,3 )\\V! -V 2 \\ W l,r. (3.13) 

ii wk wk wk ' 

Thus, if C(||?7|| L 3 + ||wi|| L 3 + II^IL 3 ) < 1, we conclude v\ = v 2 . □ 

wk wk wk 
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3.2 Pointwise bound 

In this subsection, we will prove pointwise bound of w using H^H^i,, < C*. 
For any fixed xo G Bi/ 2 \{0}, let R = \xo\/4 and Ek = B(xo, kR), k = 1,2. 
Note q* G (3, oo). Let s be the dual exponent of q*, 1/s + 1/q* — 1. We have 

IMIli^) £ lhllL^(^)l|l|l^) <C^ 4 " 3/9 - (3-14) 



By the interior estimate Lemma 12.21 

l|V^||^ (El) < WfWzsw + R- A+m M\ LHE2) (3.15) 
where / = U <g) w + w ® (U + w). Since \U\ + \w\ < C*\x\~ l < C^R' 1 in E 2 , 

\\f\\ L «* m < C^R" 1 \\w\\ Lq * (E2) < ClR~\ (3.16) 
We also have R~ 4 + 3 /i* \\w\\ Ll{Ea) < R~ 4+3 /i' 'CJ&-V q = C*R~\ Thus 

\\Vw\\ Lq * {Ei) < C*R-\ (3.17) 
By Gagliardo-Nirenberg inequality in E\, 

\\Hl°o {Ei ) ^ IhllwV) ll Vu, H^*(Si) + jR ~ 3 II^IIl 1 ^)' ( 3 - 18 ) 

where l/oo = (1 - 6>)/<f + 0(l/g* - 1/3) and thus 9 = 3/q - 1. We conclude 
— C*R~ e . Since xo is arbitrary, we have proved the pointwise bound, 
and completed the proof of Theorem 11.11 

Remark. Equivalently, one can define v(x) = u(xq + Rx), find the equation of v, 
estimate v in L°°(Bi), and then derive the bound for w(xq)- 



3.3 Proof of Theorem IT21 

In this subsection we prove Theorem ll.2l For any x G £> 2 \{0}, let v(x) = Xu(Xx+x ) 
with A = min(0.1, |x |)/2. By our choice of A, v is a very weak solution in B 2 and 
\\ v \\l 3 ,(b 2 ) — e = \\ u \\l 3 ,(b 2 i\{o})- Ml' we nave IklLoo^) — for some constant 
C 2 if £ is sufficiently small. Thus |^(x )| < C^A" 1 < 40C 2 £:|xo|~ 1 . 



3.4 Proof of Corollary PT31 

In this subsection we prove Corollary II. 5[ Suppose w satisfies (11.41) with = C*(g), 
q = 2, given in Theorem II. 1L Let 6 be given by ( 12. 8ft . U = U b and w = u — U. 
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Let u\ = Xu(Xx) be the rescaled solution and w\(x) = Xw(Xx). Note U is scaling- 

invariant. Then u\ = U + w\ also satisfies (11.41) with same C*. By Theorem 1 1 . 1 1 wit h 
q = 2, we have the bound 

\w x (x)\ < CCJx\~ 1/2 , |x| < 1, (3.19) 
which is uniform in A. In terms of w and y = Ax, we get 

\w(y)\ < CCX^X^y^ 1 ' 2 , \y\ < X. (3.20) 
Now fix y and let A — > oo. We conclude w = 0. 
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